Metric deformation and eigenvalue problem in 2D for an irregular boundary 
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We prescribe a new formulation for solving the Helmholtz equation in 2D with irregular boundary. 
A suitable diffeomorphism is used to annul the asymmetries in the boundary by mapping it into 
an equivalent circle. This results in a modification of the metric in the interior of the region and 
manifests itself in the appearance of new gauge dependent source terms in the original homogeneous 
equation. The modified equation is then solved perturbatively. This method allows us to retain 
the simple form of the boundary condition at the cost of complicating the form of the original 
equation. It is seen to work reasonably well even for boundaries with large deviations from a circle 
by comparing our results with the exactly /numerically obtained ones. The Fourier representation 
of the boundary ensures the convergence of the perturbation series. 
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Introduction. — The two dimensional Helmholtz equa- 
tion appears in a wide range of physical and engineering 
problems across diverse fields — like the study of vibra- 
tion, acoustic and electromagnetic wave propagation and 
quantum mechanics. In a large class of these problems 
one is required to determine the eigenspectrum of the 
Helmholtz operator for various boundary conditions and 
geometries. Analytic closed-form solution to the eigen- 
value problem [1,] can however be obtained only for a 
restricted class of boundaries. Invoking the geometry of 
the problem by a suitable choice of co-ordinates often aids 
in finding the solutions (e.g. elliptic boundary where, 
separation of variables lead to a solution in the form of 
Mathieu functions). However, one quickly exhausts the 
list of such problems where simplification by virtue of 
using a specific co-ordinate system is possible. In most 
physical problems, one encounters boundaries which are 
far removed from such idealization like the case of the 
quantum dot. The dots are believed to be circular but 
in practice that can hardly be guaranteed. In such a sce- 
nario it is natural to consider the confining region to be 
a supercircle [2|. Further, the problem gets analytically 
intractable for irregular boundaries. 

The problem of solving Helmholtz equation for an ir- 
regular boundary has mostly been tackled using numeri- 
cal methods Q • The analytic approach towards this has 
mainly revolved around various approximation methods. 
Of these, the perturbative techniques stand out as be- 
ing the most widely used 0, 0] . The work by [j| study 
irregular domains in a general formalism using Fourier 
representation of the boundary asymmetry treated as a 
perturbation around an equivalent circle. 

We have explored an alternative approach towards 
solving the eigenvalue problem for the two dimensional 



Helmholtz operator in the interior of a region bounded 
by an irregular closed curve. The general problem is 
mapped into an equivalent problem where the bound- 
ary is a regular closed curve (for which the Helmholtz 
equation is exactly solvable) whereas the equation itself 
gets modified owing to the deformation of the metric in 
the interior. The modified equation, we see, can be writ- 
ten as the original Helmholtz equation with additional 
terms arising from perturbation in the metric. The extra 
pieces can now be treated as a perturbation to the orig- 
inal Helmholtz operator. The equation is thereby solved 
using the Schrodinger perturbation technique [6]. The 
corrections to the eigenvalues and the eigenfunctions are 
expressed in a closed form at each order of perturba- 
tion. In this approach towards solving the equation, the 
boundary conditions are specified on some known regu- 
lar curve and maintain the same simple form through- 
out. This bypasses the issue of imposing constraints on 
a boundary having a complicated geometry. 

We expect this perturbative scheme to effectively solve 
the eigenvalue problem for boundaries which reflect slight 
departure from known regular curves. We have verified 
our method against the numerically obtained solutions 
for a square (exact), an ellipse and a supercircle Q. In 
this analysis we use Fourier representation of the pertur- 
bations. This allows us to apply the method to a general 
class of continuous or discontinuous asymmetries. 

Formulation. — The homogeneous Helmholtz equation 
on a 2 dimensional flat simply connected surface S reads, 



(</ J V,V, + A: 2 ) V-> = 0, 



(1) 
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where V represents a covariant derivative. We look for 
solutions in the interior of the bounded region with the 
Dirichlet condition ip = on the boundary dS. The pa- 
rameter k 2 may be identified as fc 2 = E for the quantum 
mechanical problem of a particle confined in the region 
having energy E. 

It is convenient to work in polar coordinate system 
(r, 0) , where any closed curve satisfies the periodicity 
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condition r(9) = r{6 + 2vr). We consider a general ir- 
regular boundary of the form r = r(9). In this analysis 
we assume that the arbitrary boundary can be expressed 
as a perturbation around an effective circle (the analy- 
sis can in principle work for a deformation around any 
simple curve for which the Helmholtz equation is exactly 
solvable). We introduce new coordinates (R, a) with the 
transformation (r, 9) — ¥ (i?, a) given by 



R + ef(R,a) & 6 = a, 



(2) 



where e is a deformation parameter. This defines a diffeo- 
morphism for the entire class of well behaved functions 
f(R,a). A suitable choice of f(R,a), shall transform 
our irregular boundary into a circle of fixed radius, say 
R = i? (= ^ Jo* r{6) M) in the R - a plane. The de- 
formation of the irregular boundary to a circle changes 
the components of the underlying metric gij (r, 9) in the 
interior (g i3 (r,9) g l3 (R,a)). 

The flat background metric in the (r, 9) system is given 
by gij = diag(l,r 2 ). Under the coordinate transforma- 
tions ((2]) this takes the form 



(l + £ / (1 < 0) ) 
e /(M) (1 + e /(i,o)) 



e /(o,i) ( 1 + e /(i,o)) 
(R + ef) 2 - 



2/(o,i) 2 



where = and /t * 1 ) = §£. Here / and all 

its derivatives are expressed as functions of (R, a) (de- 
pendence on the arguments are not shown explicitly for 
brevity). We note that except for T a RR — 0, all the 
components of connection T are non- vanishing. The dif- 
feomorphism ([2]) does not induce any spurious curvature 
in the manifold (Riemann tensor, R l j k i = 0). 

The Eq. ([T]), where %j) — -0(r, 9), transforms under the 
map (r, 9) — > (R, a) to 



2 _f (0,1)^,(1,1) 



(R + ef) 2 ~ (fl + e /)2( e /(i,o) + 1 ) 

2/(o,i)^(i,o) + e /)/(i.i) _ /(o,i) ( e /(i,o) + 

' (i? + e/) 3 (e/ (1 ^ 0) + l) 2 

(i? + e/) 2 + e 2 / ( °' 1)2 l </> (1 ' 0) / (2 ' 0) e 



(i? + e/) 2 (e/(i^) + l) J 
(R + e/) 2 - e/(°' 2 ) (R + ef) + 2e 2 f^ 2 ] ^ 



+ 



(i? + e/)3(e/(M)+l) 
(i? + e/) 2 + e 2 /(°' 1 ) 2 l ^ 2 -°) 



{R + efY{ef^) + iy 



(3) 



where ^) = 0, ^) = and ^ = 
Here ip and all its derivatives are also expressed as func- 
tions of (R, a). 

The analysis can proceed from here for a specific form 
of the function f(R,a). We choose f(R,a) = Rg(a), 
where g(a) can be expanded without a loss of generality 



in a Fourier series. We further impose g(a) = g(—a) for 
simplicity whereby only the cosine terms are retained 



g(a) = C„ cos 



na. 



(4) 



n=l 



The constant part Co can always be absorbed in R ([2]). 
With this choice of f(R,a), Eq. ^ simplifies to 



J2 * n £ni> + Ei> = 0, 



(5) 



where the operator C n is given by 



6i? 2 y 
+ n(n - l)Rg' 2 |2^ (1 

+ 6. 9 2 {^°^+i?V (1 ' 0) 



3ni? 5 { 3 , V( 1 '°)+2« ? V( 1 ^ 1 )} 



o) 



1? 



} 

ty(2,0)j 



(6) 



In the above equation a prime denotes a derivative with 
respect to a. 

We shall adopt the method of stationary perturbation 
theory [6| to solve for tp and E. Thereby, treating e as a 
perturbation parameter we expand the eigenfunction ip 
corresponding to the eigenvalue E as 



i, = V (0) + «A (1) + eV 2) + ■ • • 
E = E^+eE^+e 2 E^ + --- 



(7) 



with superscripts denoting the order of perturbation. 
We assume that the perturbative scheme converges and 
(ip,E) can be calculated order by order up to any arbi- 
trarily desired precision. We note that the parameter e 
is arbitrarily invoked to track different orders and could 
be absorbed in the Fourier coefficients C n . 

Plugging ([7]) in (O, and collecting the coefficients for 
different powers of e yields 

O(e ): (£o + £ (0) )V; (0) =0 (8a) 
©(e 1 ): (£ + £ (0) )^ (1) + (£i+£ (1) )^ (0) =0 (8b) 



C(e m ) : J2 ( £ « + E(n) ) ^ (m ~ n) 

n=0 



(8c) 



The change in the metric components induced by 
the smooth deformation @ amounts to a gauge trans- 
formation and generates source terms to the unper- 
turbed homogeneous Helmholtz equation at each or- 
der in e. At the i th order we have the terms 
Ci^°\ Citp^ ■ ■ ■ , CiTpd^ 1 ^ , which have no physical 
origin and are merely artifacts of the chosen gauge. Main- 
taining the simplicity of Dirichlet boundary condition is 
hence achieved at the cost of new terms appearing in the 
original equation. 
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The unperturbed energy E^> and corrections 
are given by 

E (0) = _ (^(o)|£ |^(o)). E (D = _^(0)| A |^(0) ); 

= - (V (0) |A + - (V> (0) |4# (0) >. (9) 

A unique feature of our method is that the Dirichlct 
boundary condition, ipw (Rq,oi) = maintains its sim- 
ple form separately for every order in perturbation. 
The general solution of the Eq. ([8a) is 



1>i°J = No*JoQ>); 
= Ni tj Ji(p) 



cos(Za) 
sin(Za) 



(1 = 0) 
(1^0) 



(10) 



where Ji is the I order Bessel function with the ar- 
gument p = ij E^ ) R and Nij is a suitable normali- 
sation constant with I g N, j G N>o- The energy 
spectrum -Ej^ of the unperturbed Helmholtz equation 
is dictated by the the j th zero of J; denoted by p l . , 

Using ^(R ,a) = 0, we have = P^/Rq- In 

the order of ascending magnitude the energy levels are 

E o°l ' ^i°i ' E 2°l i -^o°2 i ^34 ) i E fl " " ) where all 
the levels with non-zero Z are doubly degenerate. 

In this formulation the energy corrections can be ob- 

(i) 

tained in two ways. Firstly, E\ '■ can be estimated from 

the knowledge of ip^iy m < i) using Eqs. ([9]). Alter- 
natively it can be extracted by imposing the boundary 
condition, ipfj(Ro,a) = which in addition yields the 

coefficients of Bessel function (in ip^. ) . The method can 
in principle be used to calculate corrections at all orders 
of perturbation. We calculate the energy corrections for 
the following two cases. 

Case I: N on- degenerate states (I = 0) 
The first order correction to the eigenfunction is obtained 
from the Eq. (I8bl) . Thus, we have 

E iX) 

1 °4 n °Mp) 



ipoj =a Jo(p) 



2E, 



+ 



E 

P =i 



C P Ji(p)} cos (pa) ( 



11) 



The first order energy corrections are obtained by impos- 
ing the boundary condition ip^j (R Q , a) = 0. This yields 

E o] = and a p = PojNojCpJiip^yjpfaj), (p 7^0). 
The vanishing of the first order correction is verified using 
Eq. The remaining constant a can be determined 
by normalising the corrected wavef unction. 

The first non- vanishing energy correction occurs at the 

(2) 

second order. The correction Eq • is obtained by substi- 
tuting ipQj and ipQj in Eq. ^ . We have 

Po, 3 Jn(Po.j) 



P (2) _ F (o) t r 2. t ,_; , 



Jn(p 0lj 



(12) 



We may as well solve (I8cl) for m — 2 to obtain 



oE 

^ 0,7 11 — 1 



p 7 n— 1 
oo 

+ E 

P =i 



) Jn, 3 {p) cos(pa) 



^b p J p (p) - pa C p Ji(p)} cos(pa) (13) 



where J nJ (p) = §{a n J' n (p) - § N 0tj C n J{ (p)}. 

Boundary condition -0g j(i? o ,a) = extracts EqJ as 
given in (Eq. (|12p) and the coefficients 6 p . 

Case JJ; Degenerate states (I ^ 0) 
In this case, the first order wavefunction correction is 
given by 



(i) 



(I) 1 
1,3 



E, 



(0) 



cos(Za) 



|a p J p (p) + ^N h jJ[(p) (Ci+p + q ; _ p! )| cos(pa) 



aoJb(p) 



■//(p) (14) 



Here we have considered only the 'cosine' form of ipj°) 
(see Eq. (TTT))) ') for the I ^ case. The other solution with 
the sin(Za) term can be treated similarly. We estimate 
the first order energy corrections by imposing the bound- 
ary condition on ipfj . E^y* is also verified by substituting 
,,(o) 



ip^J in Eq. ©. We have 



E 



(i) _ p (o) r 



(15) 



which is generally non-vanishing unlike the earlier case. 
The second order energy correction becomes crucial when 
C21 = 0. However, in this work we have not considered 
this. The choice of the 'sine' solution for ip^ (Eq. ITOj) 



17.(1) 

gives E\J 



The coefficients ao and a p are obtained by using the 



boundary condition on ipj . '. This gives us 
N ltj p l:j d J>( Pl J 



a 



2 Mp^Y 



ZJpiPiJ 



[c p+l + qp_,|] , (p?o,i). 



The coefficient ai can be obtained from the normalisa- 
tion of ipij- These results are consistent with the earlier 
investigations @. 

Results and Discussions. — We next apply the formal- 
ism developed in the earlier section to a few specific 
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Energy labels 




Square 






Ellipse 




Supercircle 








a 




r = 


a 






a 






r 


(| COS S| + | £ 


infl|) 




Jl-E 2 COS 2 e 


r = 71 

( | cos 


9\* + \ sin0|* 








a = 1 




a = 73/2 e 


= 0.5 


a — 


1 t = 3 






H: = 0.7935 


Ro = 0.9294 


Ro ■■ 


= 1.0596 




l,k 


Es 


Ps 


% Error 


Ns 


Ps 


% Error 


Ns 


Ps % Error 


0, 1 


9.870 


10.129 


2.61 


6.744 


6.744 


0.00 


5.219 


5.217 


0.04 


U, 1 


49.348 


49.533 


0.37 


37.191 


37.777 


1.58 


27.157 


27.138 


0.07 


2, 1 


39.478 


36.047 


8.69 


29.080 


30.416 


4.59 


22.372 


22.141 


1.03 


2, 1 


49.348 


47.726 


3.29 


30.725 


30.652 


0.23 


25.087 


24.838 


0.99 


4, 1 


88.826 


87.429 


1.57 


65.637 


66.664 


1.56 


51.782 


51.365 


0.81 


4, 1 


98.696 


95.471 


3.27 


65.343 


66.662 


2.02 


50.967 


51.203 


0.46 



Table I: Comparison of the first few energy eigenvalues (with the magnitude of % error) for different boundary geometries 



boundary geometries. The analysis is restricted to con- 
vex domains only to ensure numerical convergence. We 
have considered the case of a square, an ellipse and a 
supercircle. The specific form of r(9) for these closed 
curves are used to calculate the metric deformation and 
thereby estimate energy corrections. The perturbative 
prescription is seen to converge with dominant non-zero 
corrections coming from the first few orders. It can be 
seen that the first and second order corrections in energy 
are linear and bi-linear in C n respectively. Further, it is 
clear that E^ m ' will be to- linear in C n , hence the con- 
vergence of the Fourier coefficients, C n , will ensure the 
convergence of the series. In this analytic formalism the 
approximation appears only through truncation of the 
series given by Eqs. (J7|. 

We have estimated energy corrections up to the sec- 
ond order in perturbation for the I = states and only 
first order corrections for the I ^ states. The two- 
fold degeneracy of the original I ^ states split at the 
first order for C21 7^ 0. Our results for first few energy 
levels are summarised in the Table |U The results ob- 
tained using the perturbative scheme (Ps) are compared 
with the numerical (Ns) or exact (Es) solutions. Numer- 



ical solutions are obtained using the PDE Toolbox of 
MATLAB®. Our results are in good agreement with the 
numerical results. The discrepancy is ~ 5% for the case 
of the square and < 1% for the supercircle. The method 
is expected to yield better results for smooth boundaries 
without vertices. The results for the supercircle even 
at the first order show better accuracy than its second 
order counterpart for the square. However, our method 
can achieve any desired accuracy by incorporating higher 
order corrections. 

In conclusion we note that Fourier decomposition of 
the boundary asymmetry makes the method completely 
general and holds good for a wide variety of boundaries. 
Since our solutions are general in principle other bound- 
ary conditions such as, Neumann, could also be applied 
to obtain the corresponding spectrum. This makes it 
crucially important in estimating energy eigenvalues and 
eigenfunctions for problems. 
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